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We discuss QCD with two light flavors at large baryon chemical potential /i. Color superconductivity 
leads to partial breaking of the color SU(3) group. We show that the infrared physics is governed 
by the gluodynamics of the remaining SU(2) group with an exponentially soft confinement scale 
^QCD ~ ^ exp(— a/i/(5rA)), where A ^ /i is the superconducting gap, g is the strong coupling, and 
a = 2\/2-k/11. We estimate that, at moderate baryon densities, Aqqq is 0(10 MeV) or smaller. 
The confinement radius increases exponentially with density, leading to "asymptotic deconfinement." 
The velocity of the SU(2) gluons is small due to the large dielectric constant of the medium. 



Introduction. — Soon after the discovery of asymptotic 
freedom in QCD |^ a hypothesis was put forward that, 
at high baryon densities, quarks (which are normally con- 
fined in hadrons by strong forces) are liberated, i.e., nu- 
clear matter transforms into deconfined quark matter . 
In recent years, our knowledge of dense quark matter has 
considerably expanded. We now understand that, in re- 
ality, dense matter shows more intricate features than in 
the original picture of In particular, quark matter 
at high densities exhibits the phenomenon of color su- 
perconductivity , which determines the symmetry of 
the ground state and the infrared dynamics. 

The number of light quark flavors Nf turns out to play 
a crucial role. The simplest case is Nf = 2, where up and 
down quarks are massless and other quarks are neglected. 
The following picture emerges in perturbation theory, as 
well as in instanton-inspired models. The condensation of 
color antitriplet up-down diquarks breaks the color SU(3) 
down to an SU(2) subgroup. Thus, five of the original 
eight gluons acquire "masses" by the Meissner effect [^|| , 
similar to the Higgs mechanism. The remaining three 
gluons are massless (perturbatively) . Because of Cooper 
pairing, the spectrum of quark excitations carrying SU(2) 
color charge has a gap A. 

In order to understand the physics below the energy 
scale A we must examine the pure gluodynamics in the 
remaining unbroken SU(2) sector. As we shall see, the 
process of high-density "deconfinement" is quite non- 
trivial in this case: the quarks are always confined (as- 
suming that SU(2) Yang-Mills theory confines), how- 
ever, the confinement radius grows exponentially with in- 
creasing density. We shall also see that, at scales much 
shorter than the confinement radius, the dynamics of the 
SU(2) gluons is similar to electrodynamics in a dielectric 
medium with large refraction index. 

The effective Lagrangian. — Below the scale A, we ex- 
pect that the heavy (gapped) degrees of freedom decouple 
and the remaining fields can be described by a local effec- 
tive Lagrangian. The absence of quarks carrying SU(2) 



charges below A implies that the medium is transparent 
to the SU(2) gluons: there is no Debye screening and 
Meissner effect for these gluons. Mathematically, the po- 
larization tensor nj^j (g) vanishes at g = 0, which can be 
checked by a direct calculation of 11 at small q, as was 
done in Ref. . The absence of Debye screening means 
that a static color charge inserted into the medium can- 
not be completely screened as it is in hot plasmas. This is 
easy to understand since all quarks carrying SU(2) color 
are bound into SU(2) singlet Cooper pairs. Analogously, 
the Meissner effect is absent because the superconducting 
currents, which are coherent motions of the condensate, 
cannot screen the magnetic field, since the condensate is 
SU(2) neutral. Thus, at first sight, it might seem that the 
quarks in the medium have no effect on the gluon effec- 
tive Lagrangian, which must be simply L — ~ F^^ / {4:g'^) , 
i.e., the SU(2) Yang-Mills Lagrangian with the coupling 
g matching the running coupling in the original theory 
at the scale A [||. However, a closer look shows that 
the situation is somewhat more complicated and, in fact, 
more interesting. 

Although a static SU(2) charge cannot be completely 
Debye screened by SU(2) neutral Cooper pairs, it can 
still be partially screened if the medium is polarizable, 
i.e., if it has a dielectric constant e different from unity. 
If e > 1, then the Coulomb potential between two static 
color charges is g^/{er); i.e., the gauge coupling is effec- 
tively reduced by a factor of e^/^. As explained in more 
detail below, this is exactly the situation in the color- 
superconducting phase. Analogously the medium can, in 
principle, have a magnetic permeability A ^ 1. (We de- 
note the permeability by A instead of the more common 
fi, since the latter symbol is already used for the chemical 
potential.) The dynamics of gluons is thus modified by 
the dielectric constant and the magnetic permeability of 
the medium. Hence, one needs to develop the theory of 
"gluodynamics of continuous media" , which, as far as we 
know, has never been encountered before. This theory, 
in contrast to its U(l) counterpart (the electrodynamics 
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of continuous media), is an interacting theory. 

Fortunately, even without exphcit calculation we can 
already write the effective Lagrangian of SU(2) gluons 
from rather general arguments. It should satisfy the re- 
quirements of locality (since the quarks that have been 
integrated out have gaps) and gauge invariance. It does 
not need to be Lorentz invariant, since this invariance 
is already violated by the presence of the high-density 
medium, but it should be rotationally invariant and con- 
serve parity. Thus the effective action at the scale A 
must have the following form 



'S'off ~ 
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and B°- = ^(-ijkFj^.. Higher-order cor- 



where Ef = F^, 

rections (in powers of fields and derivatives) to (□) are ir- 
relevant for the infrared physics and have been neglected. 
The constants e and A in Eq. (|^) have the meaning of the 
the dielectric constant and the magnetic permeability in 
the regime where the gluon fields are linear. The speed 
of gluons in this regime is u = 1/ \/eA. 

As we shall see, eA ^ 1. This should be contrasted 
with the standard picture of the vacuum as a dielectric, 
in which case eoAo = 1 by Lorentz invariance. In other 
words, vacuum polarization effects (such as gluon, ghost, 
and /X = quark loops) in Eq. (|l|) have been absorbed 
into the running of . The constants e and A are defined 
relative to those of the vacuum. 

One can give a parametrically correct estimate of e 
using a crude model, in which the Cooper pairs are rep- 
resented by classical oscillators with the spring constant 
k. The characteristic scale for a Cooper pair "oscillator" 
is A; thus one can estimate A: '--^ A'^. The polarizabil- 
ity of each oscillator is g^ /k. Multiplying by the density 
of the Cooper pairs, /l(^A, one finds e — 1 ~ g^/i^/A^. 
This estimate is corroborated by the explicit calculation, 
presented later in this paper, where we obtain 



e = 1 
A = 1, 



K = 1 
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At high densities, the gap A is exponentially sup- 
pressed compared to the chemical potential ^ 0, 
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where g is the gauge coupling at the scale /i, and h is 
some numerical constant. According to Eq. (^), k ^ 1 
and we can write 
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which means that the dielectric constant of the medium is 
very large. Hence, the Coulomb potential between SU(2) 



color charges is greatly reduced. This can be interpreted 
as a consequence of the fact that the Cooper pairs have 
large size (of order 1/ A) and so are easy to polarize. The 
magnetic permeability, in contrast, remains close to 1 due 
to the absence of mechanisms that would strongly screen 
the magnetic field. 

The scale of confinement. — Once the effective La- 
grangian (|l|) is obtained, one can use it to investigate 
the infrared dynamics of the gluons. One notices that 
(|l|) possesses a modified Lorentz symmetry in which the 
speed of light c = 1 is replaced by 
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One can make this symmetry manifest by rescaling the 
time, the field Aq and the coupling in Eq. ([^), 



(7) 



After the rescaling (^, the action (|l|) assumes the famil- 
iar Lorentz-invariant form in the new coordinates. 



where 



F?. 



(8) 
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The coupling in the action (^ is not g but g' which is 
smaller by a factor of e^l^. This means that the small 
parameter that controls the perturbative expansion in 
the theory (||) is not as = g^ li^ii) but rather 



(10) 



which is much smaller than as, since e is large. 

Another way to derive Eq. ( p^ ) is by restoring h and 
c in the expression for the strong coupling constant a^, 
which is given by g^ j {Airhc) in the vacuum. In our dielec- 
tric medium, the Coulomb potential between two static 
charges separated by r is g^/{er). Thus, we have to re- 
place g^ by g^g = g^/e. The velocity of light c also needs 
to be replaced by the velocity of gluons v. This gives 
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(11) 



since ?i = 1 in our unit system. Equation ( [TTI ) coincides 
with in Eq. ([l0|), as one expects. 

Using Eq. d), one can express the coupling a'^ in terms 
of the gap A, 



3 gA 

*s = 1= 

2x/2 [I 



(12) 
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Equations ( p^ and ( [T^ define the couphng in our effec- 
tive theory at the matching scale with the original micro- 
scopic theory, i.e., at the scale A. The coupling increases 
logarithmically as one moves to lower energies, since pure 
SU(2) Yang-Mills theory is asymptotically free. This 
coupling becomes large at the confinement scale Aq^j-,, 
which is the mass scale of SU(2) glueballs. The spectrum 
of these glueballs is known from lattice studies of SU(2) 
Yang-Mills theory ||], except that the role of the speed 
of light is now played by v [Eq. (||)]. Since a'^ is tiny 
[because A//x <C 1 in Eq. (p^)] it takes long to grow, and 
the scale Aqq^ is thus very small. Using the one-loop 
beta function, one can estimate 
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where (3q is the first coefficient in the beta function and 
is equal to 22/3 in SU(2) gluodynamics. 

We can draw a few immediate conclusions from Eq. 
(p^). First, Aq(-,p depends very sensitively on the gap 
A, in particular on the numerical value of the constant b 
in Eq. (||). Unfortunately, the latter is not exactly known. 
The uncertainty in the value of the gap A translates into 
a huge variation of Aq^-,]-,. For example, if one uses 



b= dUn^ 



(14) 



which is obtained by solving the one-loop gap equation 
where the exchanged gluon propagator is replaced by the 
hard dense loop (HDL) expression ^ , then with Aqcd = 



200 MeV we find A' 



However, if we use 



QCD 



b = 5127r exp 



10 MeV at /I = 600 MeV. 



(15) 



which is obtained if one assumes the Bardeen-Cooper- 
Schrieffer ratio between the critical temperature Tc and 
the gap A, and computes Tc by taking into account the 
fermion wave- function renormalization ||l0[] , then, at the 
same chemical potential, Aqg^^is reduced to a mere 0.3 
keV. Regretfully, neither Eq. (|lj) nor ( |l5| ) seems to be en- 
tirely correct, since there are physical effects that they do 
not take into account (e.g., the Meissner effect). Clearly, 
any attempt to give even the roughest numerical estimate 
for Aqqj-, requires an accurate determination of the gap 
A. It has been argued that to compute A one needs a 
better understanding of the issue of gauge invariance, the 
finite fermion lifetime, and the running of the coupling 
| pj| . Regardless of all these uncertainties, the exponen- 
tial dependence of Aq(-.]3 on /i/A makes it safe to predict 
that, even at moderate values of fi, the confinement scale 
-^QCD very small, much smaller than Aqcd- 

Second, as the density, i.e., ^, is increased, A'qqq van- 
ishes exponentially fast due to the factor /i/A in the ex- 
ponent. We arrive at the following picture of how de- 
confinement occurs at large densities. Strictly speaking, 



at any given density, the theory is confined. However, 
the confinement radius 1/Aqqj-, grows exponentially as 
the density is increased. Therefore, if one looks at the 
physics at some large, but fixed, distance scale, there is 
a crossover density when effectively the color degrees of 
freedom become deconfined at that scale. We call this 
phenomenon "asymptotic deconfinement." 

The computation of e and A. — To find e and A, one 
has to calculate (|l|) by integrating out the quark degrees 
of freedom in the QCD Lagrangian. This amounts to 
computing the one-loop polarization operator Il{q) and 
the gluon vertices r3(gi, 92), ^i{qi, 92, qs), etc. This pro- 
cedure is the same as the one giving rise to the hard 
thermal loop (HTL) and hard dense loop effective ac- 
tions ||l2| . The situation here is simpler than in the HTL 
and HDL cases: in the regime where all gluon momenta 
q are much smaller than A, the functions H and F can 
be expanded in powers of q, yielding a local effective La- 
grangian. (In contrast, the HTL and HDL actions are 
non-local, since the fermions do not have gaps.) The 
gauge invariance of the effective Lagrangian greatly sim- 
plifies our task: in order to know e and A, one needs 
to compute only the polarization tensor H of the SU(2) 
gluons. The leading contribution at large density comes 
from the superconducting quark loop. The detailed cal- 
culation of H was done in Rcf. From Eq. (99a) of 
that paper one can derive the following expression for 
nOO(Qo,q),a,6=l,2,3: 

00 q/2A 

n°g(qo,q) = -5a.^- fdz f dy(l ' + ^ 
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u+ + u-+go/A + — go/A 



(16) 



where u± — •\/(z ± y)"^ -\- 1, and we assume that go a-nd 
q = |q| are much smaller than so that the dominant 
contribution comes only from particles near the Fermi 
surface. Physically, the first term in parentheses is the 
probability to excite a quark-hole pair through an SU(2) 
gluon, and the last term contains the corresponding en- 
ergy denominators for such an excitation. For jgol > 2A, 
one should replace go by qo + ie. Expanding Eq. ( p^ ) to 
quadratic order in qo and q around qo = q = 0, one finds 



H 



where 



afc('?0,q) = -Hq Sab 
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The appearance of a factor /i^ is due to the fact that 
the loop integral is dominated by the momentum region 
near the Fermi surface, whose area is proportional to /i^. 
Similarly, one obtains from Eq. (99b) of ||] the following: 



n°U90,q) 



-nq^q^'Sab ■ 



(19) 
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The computation of n*"'jj(go, q) is faciHtated by writing 

n(go, q) = [n(qo, q) - n(0, q)] + [n(0, q) - HnDLlO, q)] 
+ nHDL(0,q), (20) 

where IIhdl is the standard HDL gluon self-energy, 
which vanishes for qo — 0. The term 11(0, q) — IIhdl (0, q) 
can be shown to be of order ©(A^) and thus negligible 
compared to the first term Il(qo, q) — 11(0, q). The reason 
for the manipulation ( po| ) is to get rid of the antiparti- 
cle contributions in this term. With Eq. (99c) of 
result for II^^j, can be written analogously to Eq. 

oo g/2A 2tt 

nrb(9o,q) = -<5a.^^ JdzJ dyj^yy 
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where k = (sin cos sin 6* sin (/s, cos 6*), and cos 6* = 
2yA/q. Expanding to quadratic order in go a-nd q, one 
finds 



n'^,{qo,Cl) = ~^qlS''Sab. 



(22) 



Note that (i) the polarization tensor satisfies current con- 
servation, g^jll^^ — 0; (ii) at this order, there is no spatial 
transverse contribution ~ q'^6^^ — q^q^ to IT^ , although 
such a term is not forbidden by the symmetries. 

After the quark loop has been integrated out, the 
quadratic term in the effective Lagrangian becomes 

An-q)[D;:H<l) + n^.(9M"('Z)/(V), where D is the 
bare gluon propagator. Comparing with the quadratic 
terms in Eq. (|^), we obtain Eqs. (||) and The fact 
that A = 1 [Eq. (H)] is due to the absence of the spatial 
transverse term in II*-' . 

In our analysis above we have neglected other low en- 
ergy excitations: the unpaired fermions of the third color 
and the pseudoscalar isoscalar mode similar to the ry me- 
son. This is justified because they are colorless with re- 
spect to the unbroken SU(2)c gluons. It is also interest- 
ing to note that since the gap in the SU(2) colored quark 
spectrum. A, is much larger than Aqqq, the spectrum 
of mesons made of these quarks must resemble that of 
heavy quarkonia. 

The asymptotic deconfinement phenomenon is spe- 
cific to the case Nf — 2. At sufficiently high densities 
when, effectively, Nf — 3, the ground state is the color- 
flavor locked (CFL) state: the color symmetry is broken 
completely [Q and all gluons are screened, both elec- 
trically and magnetically [ p^ . The low energy modes 
in the CFL phase are the Goldstone modes described 
by a chiral effective Lagrangian [|l5|,|6|. Between the 
two regimes, when one has strange quark matter below 



the "unlocking" phase transition [Q, the system might 
be an anisotropic dielectric if an ss condensate breaks 
rotational invariance. Another interesting regime with 
asymptotic deconfinement is that of high isospin density 
p8[ , where the gauge group of the gluodynamics of con- 
tinuous media is SU(3). 
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